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Abstract 



We demonstrate the close similarity of a generalized Fermi breakup model, in which densities of 
excited states are taken into account, to the microcanonical statistical multifragmentation model 



used to describe the desintegration of highly excited fragments of nuclear reactions. 
PACS numbers: 25.70.Pq, 24.60.k 
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I. INTRODUCTION 



Both the Fermi breakup (FBM) and the statistical multifragmentation (SMM) models 
provide prescriptions for calculating mass and charge distributions and multiplicities of the 
fragments emitted in the breakup of an excited nuclear system. Yet they would seem to be 
very different models. They are formulated in diferent terms and usually applied in very 
different regions of mass and excitation energy. We will show that they are much more 
closely related than they might appear to be at first glance. 

The FBM was originally proposed as a means of calculating the multiplicities^ and an- 
gular distributions^] of pions and antiprotons produced in high-energy collisions of cosmic- 
ray protons with nucleons in the atmosphere. It was found to be quite successful in this 
respect j^]. It was later applied to the statistical decay of light fragments of proton-induced 
spallation reactionsj4-6| and is now included as the preferred option for the equilibrium 
statistical decay of light fragments in widely-used nuclear reaction/transport codes, such as 



FLUKApl and GEANT4 



A variant called the phase space model, which partially takes 



incomplete equilibration into account, plays an important role in the analysis of experimental 
multi-particle fragmentation spectra in light-ion reactions {9- 13 1. In the context of nuclear 
reactions, the FBM usually assumes that fragments are emitted in their ground states or in 
(almost) particle-stable excited states and is formulated directly in terms of a phase-space 
integral limited only by the constraints of linear momentum and energy conservation. 

The SMM is used to describe the decay of highly-excited fragments of heavy-ion or 
spallation reactions. It assumes thermal equilibrium and thus allows for the emission of 



particle-unstable excited fragments consistent with that equilibrium. It has 
compared to experimental data and found to reproduce them reasonably^ we 
though many versions of the SMM have been proposed over the years 



aeen widely 



developed systematically in Refs. 



18-20 



18 



14 



17j. Al- 



26] . it was first 



The SMM is normally formulated in terms of a 
statistical partition function, be it microcanonical, canonical or grand canonical. The most 
appropriate of these is the microcanonical partition function, for which charge, mass number 
and energy are strictly conserved. The canonical and grand canonical partition functions 
are useful for deriving analytical or semi-analytical expressions that would be impossible 
to obtain in the microcanonical formulation or, when fluctuations are small, to simplify 
calculations. 
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In the following we will demonstrate how a model very similar to the SMM can be obtained 
from an appropriately generalized FBM. We will then discuss the differences between the 
SMM and the generalized FBM derived here and consider possible directions of future work. 



II. A GENERALIZED FBM 

A justification of the FBM would begin with the transition rate 7o-s. n from a state to 
the states of a configuration n, 

70-m = -jT P0-m| W n , (1) 

where uo n is the density of (linear momentum and energy conserving) states of the configura- 
tion per unit energy and, most importantly, the transition matrix element r ^ n is assumed 
to be independent of the individual momenta of the final configuration. In the steady state, 
the probability of producing a given configuration n can then be calculated as 

Pn= ^ n , (2) 

where the sum in the denominator runs over all possible configurations. When the transition 
matrix element is configuration independent as well, 

r 0^n = T i (3) 

then the probability of a configuration depends only on its density of final states, 

(4) 

0| w m Z^m^m 

This is assumed to be the case in the FBM. 

In applications of the FBM to nuclear decay j3, {J, the phase-space integral that determines 
the density of final states of a configuration of n fragments is usually written as 

k / t r \ n—1 n 
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where the sums and products j = 1, . . . ,n run over all fragments of the breakup mode, 
while the sum I = 1, ■ • • , k runs over the distinct fragments and takes into account their 



multiplicities. Here, is the excitation energy of the decaying nucleus, Bq its binding 
energy and E c q is a term associated with the Wigner-Seitz correction to the Coulomb energy 
of the system. V n is the volume in which the momentum states are normalized and is usually 
defined as js| 

V n = (l + X )V , (6) 

where Vo is the ground state volume of the decaying nucleus and the expansion factor \ is 
usually taken to be % = 1. For the fragments, Bj is the binding energy of fragment j and 
gj is its spin multiplicity, while the E C j represent the remaining Wigner-Seitz corrections to 
the Coulomb energy, taken to be 
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p _ CqouI Z j 

°'"(i + x) 1/3 Af (7) 

Conservation of nucleon number and charge requires that 

n k n k 

A = J2 A 3 = J2 NlAl and Z = ^ t Z j = ^ t N l Zi, (8) 

j=l 1=1 j=l 1=1 

where Zj and Aj are the charge and mass number, respectively, of fragment j. The FBM 
assumes that the fragments are emitted in their ground states or in (almost) particle-stable 
excited states. 

As the total excitation energy is increased, other particle-unstable excited states that 
are long-lived in comparison to the initial decaying nucleus could also be included and can 



make significant contributions to the phase space integral 27]. These can be incorporated 
compactly using the densities of excited states of the fragments. Such an extension of the 
Fermi breakup integral takes the form 

■* = n^(^)""7n^teft) <»> 



X 



f[ ( W J fe) de i) 6 y £ o -B -E c0 -J2 + e 3 ~ B i ~ E co^j 

where Bj is the excitation energy of fragment j, Uj (sj) its density of states and Bj is now its 
ground-state binding energy. Note that this expression does not contain the fragment spin 
multiplicities, gj, which are now assumed to be incorporated in the density of states. For a 
particle with no excited states, we have u>j (e 3 -) = gjS (ej). 



4 



After rewriting the densities of fragment states in terms of the internal Helmholtz free 



energies, defined for fragment j by 



25fl 



de 3 e h'iuj {ej) , 



(10) 



dp exp [-13 {F n {(5) - E )] 



(11) 



all but one of the integrals in Eq. can be performed analytically, as is shown in the 
Appendix. We can then write the density of final states u n as 

j /"C+ioo 

2m 
where 

E = e - B (12) 
and the total Helmholtz free energy F n (j3) has been defined as 

k 

F n ((3) = Y, N i {ft (P) + fi rans 09) ~Bi- E d ) - {ft ans (P) - E c0 ) , (13) 
i=i 

with the sum over fragments replaced by a sum over distinct fragments times their multi- 
plicities. The translational Helmholtz free energies are given by 



fr u iP) = - 



In Vn 



( m N Aj \ 



3/2 N 



In (JV,!) 



(14) 



\27rh 2 P J ) Ni 

where we write the mas of fragment I as m; = rrij^Ai and the mass of the decaying nucleus 
as mo = m^Ao, with the nucleon mass. We emphasize that we have made no approxi- 
mations up to this point. The expression given in Eq. (fTTl) is exactly equivalent to that of 
Eq. ©. 

To approximate the final integral, we use the method of steepest descent. We look for a 
value O for which 



±U>f.UD) 



E n = 0. 



Using the relations of the Helmholtz free energy to the entropy and energy, 



df 



df 



and e = f + Ts = f + py-, 

dp 



(15) 



(16) 



dT ' dp 

respectively, we find the saddle point condition to be equivalent to the requirement that 
energy is conserved, 



Ni (e; (P ) + ef ans (P ) - B l - E d ) - (e% ans (P ) - E c0 ) = e - B . 



(17) 



i=i 



At the saddle point /3q, the argument of the exponential is then the total entropy, S n (Po) 

k 



P (F n (P )-E ) = plY^NAjL + 



1=1 



df* dtf rans 



trans 



jftrai 
2 "id 



dp dp J u dp 
N t (a? (A)) + sj rans (A))) - 4 ms (A>) = S n (P ) 



(18) 



i=i 



To complete the evaluation, we must calculate the second derivative to determine the direc- 
tion of steepest descent. Since 



S n = —dF n /dT and dS n /dT = Cy ;n /T, 



where Cy, n is the specific heat of the configuration at constant volume, we have 

and — T (,'K (.!)) I = -Cv n T 2 . 



Cv " ' J2 (PF n (P))=T 3d2Fn 



dT 2 T dp 2 

Near the saddle point, we then find, with T = l/Po 



dT 2 



1 d 2 



P (F n (P) - E Q ) ~ -Po (F n (p Q ) -E )-- — (PF n (P)) 
« S(T ) + ^C v , n T*(P-P ) 2 . 



(P-Po 



(19) 



(20) 



(21) 



00 



We thus conclude that the direction of steepest descent is purely imaginary. The integral 
then yields 



exp (S n (T )) 



(22) 



A/27rCy iri T 2 

which is valid as long as the specific heat Cy. n is positive. When one uses the usual SMM 



approximation to the internal energy 
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23 



25] to evaluate this expression, the 



specific heat can become negative when the negative surface term of the internal energy of 
one of the fragments dominates the positive bulk term. But in that case, the fragment can 
no longer be considered to exist. 

As To 0, we expect the phase space integral to reduce to the original form of the FBM, 
Eq.([5l), for which a well-known closed-form expression exists, 



1=1 ,n j=l 



3/2 



n-1 



E 



3(n-l)/2-l 
kin 



(2vr) 3/2 h*J r (3 (n - 1) /2) 



where 



E, 



kin 



£o — B — E c0 + {Bj + E, 



<j ■ 



as T -> , (23) 



(24) 
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is the total kinetic energy of the fragments. The steepest-descent approximation does indeed 
reduce to this closed-form expression, with all fragments in their ground states, as To — > 0, 
except for a multiplicative factor R (3 (n — 1) /2) that substitutes a Stirling approximation 
for the gamma function in the denominator, 

R(h) = T (h) I exp \(h - 1/2) In (h) - h + In (2tt) /2] . 

This factor is approximately 1.06 for n = 2, 1.03 for n = 3 and decreases to one as n 
increases. 



III. COMPARISON OF THE FBM AND THE GENERALIZED FBM 



To get an idea of the differences to be expected between a calculation using the usual 
FBM and its generalized version, we compare calculations using the two models here. To 
set a context, we consider the case of a 62 MeV proton incident on 16 O, for which double 



differential emission data exist 



28|. We describe the initial stage of the reaction using 



the Monte Carlo exciton cascade model of Blann and Chadwick |29|, [30|, which provides a 
good description of the double differential data. This model simulates the pre-equilibrium 
stage of a nucleon- nucleus collison by following a cascade of particle-particle and particle- 
hole interactions together with particle emission, until all remaining nucleons have energy 
smaller than their separation energy from the residual nucleus. When it is applied to the 
case of a 62 MeV proton incident on 16 O, it furnishes a 16 O primary compound nucleus 
population corresponding to 66% of the reaction cross section of approximately 400 mb, 
with the primary populations of the 15 and 15 N compound nuclei accounting for another 
32% of the reaction cross section. The distribution in excitation energy of the compound 
populations is extremely broad, with that of 16 O extending from zero to the center-of-mass 
energy of 58.4 MeV with a peak at about 50 MeV. 

Using this as motivation, we compare the results for the Fermi breakup of 16 at an 
excitation energy of 50 MeV including 1) only the ground states of the fragments, 2) the 
ground states and particle-bound states of the fragments and 3) the ground states and all 
excited states of the fragments found in the RIPL-2 nuclear level library |3l| . We note that 
the calculation including ground states and particle-bound states is not entirely consistent 
as it includes particle-unbound ground states, such as those of 5 He and 8 Be, which make 
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Figure 1: (Color online) Primary fragment production cross sections of Be, B, C and N isotopes 
function of the neutron excess, N — Z. 



important contributions to the primary fragmentation yields. We also point out that the 
calculation using all excited states cannot be considered complete, as the discrete level sets 
of the RIPL-2 library are incomplete at the energies available to the heavier fragments of 
several two-body decay modes, in particular, of n+ 15 0, p+ 15 N and a+ 12 C. Nevertheless, 
the discrete level sets contain a sufficient number of levels to clarify the principal differences 
to be expected between the models. 

The calculations were performed using the steepest-descent approximation in all cases. 
They furnish 16 two- fragment, 38 three- fragment, 33 four-fragment, 8 five-fragment and 
2 six-fragment decay channels. The contributions of the five- and six-fragment modes are 
negligible and the contributions of the four-fragment modes are small. The primary fragment 
production cross sections of the isotopes of nitrogen, carbon, boron and beryllium are shown 
in Fig. 1 as functions of the neutron excess, N — Z. The ground-state-only calculation 
yields flatter isotopic distributions than the others, which, with the exception of boron, do 
not have their peaks at iV = Z, the minimum of the valley of stability. The particle-stable 
state calculation tends to be more irregular than the ground-state-only one. It depresses 8 Be 
production, due to its lack of particle-bound states to compete with those of other nuclei. It 
places the nitrogen isotope peak at A = 14, but maintains the carbon isotope maximum at 
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Figure 2: (Color online) Primary fragment production cross sections of H and He isotopes as a 
function of mass number. 

A = 11, as the n + a+ n C decay mode competes favorably with the a+ 12 C mode due to the 
limited number of particle-bound excited states of 12 C. When all excited states are included 
in the calculation, all isotopic distributions have their peaks at the minimum of the valley 
of stability and roughly reflect the (inverted) form of the valley. 

The mean primary fragment multiplicity decreases from 2.8 for the ground-state only 
calculation to 2.4 for the particle-stable state one to 2.3 when all excited states are included. 
Neglecting the four- fragment or higher modes, which correspond to less than 1% of the cross 
section in all cases, these multiplicities imply that the contribution of three-fragment modes 
to the primary fragment distribution is about 80% in the ground-state only calculation, about 
40% in the particle-stable calculation and only 25% in the calculation containing all excited 
states. That is, the decay tends from dominance of the three-fragment modes to dominance 
of the two-fragment modes as excited states are taken into account. The decrease in the 
primary multiplicity is easily understood, since an increasing portion of the excitation energy 
remains in the fragments as excited states are added, rather than being liberated as a larger 
number of smaller, lesser-bound fragments. The fact that the multiplicity decreases from a 
value close to three to one close to two is a result of the relatively low excitation energy of the 
calculation. At higher excitation energies, the average multiplicity is larger but follows the 
same trend - as more excited states are included, the average primary multiplicity decreases. 
We note that the excited fragments of the generalized FBM will, of course, subsequently 
decay and increase the net multiplicity. To be consistent, however, this subsequent decay 
would also be best described using the generalized version of the FBM. 

Before concluding this section, we briefly discuss the primary distributions of hydrogen 
and helium isotopes, shown in Fig. 2. These remain fairly stable in relative yield but 
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decrease slightly in magnitude as excited states are included. The reduction in magnitude is 
a reflection of the decrease in multiplicity - the lower multiplicity of the calculations including 
excited states is mainly due to decreased emission of H and He fragments. The similarity 
in relative yields of these isotopes is to be expected, as they have either no excited states 
or very few excited states at high energy that cannot increase their production greatly, but 
are emitted in fragmentations that leave other nuclei in excited states. An exception to 
this argument is seen in the production cross section of 4 H, which comes exclusively from 
the two-fragment ground-state decay mode of 4 H+ 12 N. As neither of these residual nuclei 
have known excited states, this mode is suppressed as the excited states of other nuclei are 
included. 

Two further comments should be made in reference to calculations using the generalized 
FBM. As alluded to above, sets of known discrete levels are limited, even for stable nuclei. 
Continuum level densities are thus needed to perform reasonably realistic calculations at 
energies higher than those we have shown here. Level density parameters that have been 
fit to discrete levels and resonance densities can be found in the RIPL-2 library, but only 
for 20 F and heavier nuclei. The parameters for lighter nuclei must be obtained from the 
extrapolation of systematics or from theoretical calculations. Second, hot nuclei have a 
limiting temperature/excitation energy, above which they no longer exist. This should 
be taken into account in FBM calculations, especially in the case of light nuclei far from 
stability, where this temperature/excitation energy is expected to be quite low. One manner 
of obtaining estimates to both the level densities and the limiting temperatures of arbitrary 
nuclei is through the use of self-consistent temperature-dependent mean field calculations, 
such as those we have recently applied to the SMM 26|. 



IV. FROM THE FBM TO THE SMM 



The generalized FBM is very similar in form to the SMM used to describe the fragmen- 
tation of highly excited heavy nuclei. We discuss the similarities and differences of the two 
models here. To be brief, we refer to the generalized FBM as simply the FBM in this section. 
As many versions of the SMM have been proposed and we do not wish to compare the FBM 
to all of them, we will take as our reference the microcanonical SMM presented in Ref. ^25 1 . 
which we will denote simply as the SMM. 
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We begin by observing that the SMM uses the Helmholtz free energy of Eq. (TT3j) to 
define the entropy of Eq. (TPS]) and uses the energy condition of Eq. (fT7|) to determine the 
configuration temperature T . It then defines the statistical weight of a configuration as 

w n Ae = exp (S n (T )) , (25) 

where Aeq is a small interval in energy about the total value Eq. In most cases, the differences 
between the statistical weight of the FBM, Eq.©, and that of the SMM, Eq. are 
irrelevant, as the variations in specific heats and temperatures among configurations are 



25|]. 



small compared to the exponential variation with respect to the entropy 

In fact, although we entered into great detail to calculate the density of states of an 
FBM fragment configuration in terms of the quantities used in the SMM, the association 
between the two models was already established once we identified the statistical weight of 
the FBM as the density of states. The microcanonical SMM defines the statistical weight 
of a fragment configuration in terms of the microcanonical entropy S n , which, in turn, can 
be defined in terms of the density of states u n as 

S n = \n(u n Ae ) . (26) 

Thus, the FBM and the SMM use the same physical quantities in much the same way. They 
still have their differences, however, which we discuss next. 

A minor distinction between the FBM and the SMM, as they are commonly used, is 
their treatment of the normalization volume V n . Both models use a normalization volume 
larger than the volume of the decaying nucleus Vo, with the difference expressed in terms of 
a multiplicative factor x- As mentioned before, calculations with the FBM model often use 
a normalization volume twice that of the volume of the decaying nucleusjsj, 

V n , F BM = (1+x)Vq with X = 1 • (27) 
SMM calculations normally use a factor of \ between 2 and 5, but exclude the volume of 



the fragments 



18 



25] . so that 



V n ,sMM = y(l + X)V -J2 V jj = * V ° with X = 2 - 5. (28) 

As defined here, the volumes in both models are independent of the fragment configuration 
and are obviously the same when x = 2 is used in the SMM. In Ref. 19, a multiplicity- 
dependent volume was used in order to assure that the fragments were formed outside 
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their respective Coulomb radii. This effectively introduced a factor \ that increases with 
the total excitation energy, thereby increasing the volume available to the fragments as 
the excitation energy increases. Recently, we have used temperature-dependent fragment 
volumes Vj obtained from self-consistent calculations of the structure of hot nuclei in the 



SMM 



26|. In this case, the normalization volume is temperature and configuration dependent 



and no longer reduces to x^o- We note that the reduction of the normalization volume due 
to exclusion of the fragment volumes can be justified in both the FBM and the SMM, based 
on considerations similar to those used to obtain the van der Waals approximation to the 
equation of state of a near-ideal gas. Numerical calculations show the volume reduction to 



32 



33], but to be reasonably well 



depend on both the masses and the number of fragments 
described by the approximate form given above. 

The most important difference between the two models could be considered one of phi- 
losophy. The FBM is a model of nuclear decay while the SMM is an equilibrium statistical 
model whose configurations are identified with the fragmentation modes of the decaying 
nucleus. This difference is reflected in the fact that the SMM considers the configuration 
containing only one fragment, the decaying nucleus, that the FBM does not take into ac- 
count. This has been justified by characterizing the SMM decay as explosive and contrasting 
it to the slower compound nucleus (CN) decay, which all residual fragments are assumed to 
undergo, including the remaining fraction of the original (one-fragment) configuration, after 
the initial fragmentation^, 23 1. 

Unfortunately, neither the FBM nor the SMM furnish decay widths or lifetimes that 
could be used to compare their characteristic time scales with those of CN decay. However, 
one property accessible in both the FBM/SMM and the CN decay models is the average 
energy of the emitted particles. In the SMM, it has been shown that collective flow due 
to radial expansion contributes little to the fragment energies [34]. The average relative 
asymptotic energy of the fragments of a two-body FBM/SMM decay (assuming fragment 
volumes independent of the temperature) can then be taken to be 3To/2 + V c , where V c is 
the energy gained due to the post-emission Coulomb repulsion of the two fragments. The 
Weisskopf approximation to CN emission of a particle of type c (two-body decay) furnishes 
a statistical weight that can be written as, 

f°° 2u c e c 

27rp (e ) T c (e ) = J de c g c —^o c (e c ) p c (e - e c - Q c ) (29) 
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where po and p c are the level densities of the initial and residual nuclei, respectively, Q c is 
the Q- value of the reaction and cr c (e c ) is the absorption cross section for particles of type 
c incident on the residual nucleus at energy e c . This implies an average relative energy for 
emission of particles of type c of 2T C + V c , where V c is the effective Coulomb barrier between 
the emitted particle c and the residual nucleus and the temperature T c is very close to the 
temperature To obtained from the same two-body decay in the FBM [27]. Judging from 
the energy released by the two-body decay mode, we would thus have to conclude that 
the FBM/SMM emission is no more explosive than the CN one. We thus suspect that the 
distinction made between FBM/SMM decay modes and those of the CN is a spurious one. 
If this is the case, there is no reason to retain the original one-body configuration in the 
SMM, to later decay by CN emission, as this emission is already taken into account by the 
FBM/SMM two-body decay modes. 



V. CONCLUDING REMARKS 



We have shown that the FBM and the microcanonical SMM can be considered to be 
essentially one and the same model, if the FBM is generalized to include excited states 
and the one-fragment configuration is excluded from the SMM. The sequential CN decay 
assumed to occur after fragmentation in the SMM would then also be described in a more 
consistent manner by application of the FBM/SMM itself to the fragments. That is, the 
sequential CN decay would be substituted by sequential multifragmentation. This would 
naturally resolve the arbitrary division imposed in many calculations of post-fragmentation 
decay through sequential two-body CN decay for heavy nuclei and Fermi breakup for light 
nuclei. It, however, has the drawback of replacing well-known and very well-established 
expressions for two-body decay, such as the Weisskopf one of Eq. (1291) . with the simple 
FBM/SMM expression. The ideal solution would be an improved FBM/SMM, in which the 
two-body decay is described by a Weisskopf-like expression and n-body decay by an appro- 
priate extension of this. We would then have a consistent model of equilibrium statistical 
decay, in which the decay modes are dictated by the available energy and the characteristics 
of the system rather than by its modelers. Work in this direction is in progress. 
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Appendix A: Evaluation of the phase space integral 

We wish to evaluate the final density of states given by the extended Fermi breakup 
integral of Eq.Q, which is repeated below, 

k 1 / T; x n-l f n 



„ n / n , 

J II ( w i ^ de ^ 6 y £ ° - B o~ E co~Yl [~< 



x / I I ( w i ( £ j) d£ j) $\ £ o-B -E c0 -^[ + ej - Bj - E cj 



We begin by using the formal relation between the densities of states Uj (ej) and the corre- 



f oo 



sponding internal Helmholtz free energies /* ((3j) 25] 

' deje-teuj . (A2) 
We invert the Laplace transform to obtain an expression for the density of states ujj (sj) as 

u s ( £j ) = — / d^e^e~W^\ (A3) 



o 



2vri 

where Cj is a positive number to the right of all singularities in the complex plane. Since 
the density of states Uj (Ej) normally grows as exp (2y/a£j), the constant Cj in the inverse 
Laplace transform could be taken to zero without any effect on the result. (The integral 
defining /* (/3j) converges for any (3j with 3ft [/3j] > 0.) However, to facilitate the evaluation 
of intermediate results, it is convenient to leave it free for the moment. Substituting for the 
densities, we have 

k 1 n / \r \ n ~ 1 r n ( n \ /l\ nra / rci+ioo 

- - n^n^y /n^,(p) (^ n(/ 
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(A4) 



We begin by integrating over the excitation energies. The first integral, over e±, for example, 
furnishes 



n^n(^) /n-Us* fe)sU 
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x exp 



i=i 



— - — — Ba — E r i 

2 mj 3 CJ 



n POO 

U / e^-^dej .(A5) 

j=2 ^° 



The remaining integrals over e 2> - ■ - ,s n can now be performed if we take 3ft [fa] > 9ft [fa], 
j = 2, ... ,n. The result is 



n^n 
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(A6) 



i - Pj 



We next perform the integrals over fa, j — 2, . . . ,n. To do this, we use the fact that the 
function exp [—faff (fa)] is well behaved (no singularities) for 9ft [fa] > 0. In particular, we 
note that as 9ft [fa] — > oo, exp [—faff (fa)] — > 9j, where gj is the ground-state degeneracy of 
fragment j. We can thus close the fa contour to the right, to obtain 



1 

2ni 



Cj +ioc e -tof;<fii) 



dfa 



(A7) 



fa ~ fa 

since 9ft [fa] > 9ft [fa]. After performing the fa integrals, j = 2, ...,n, we have (taking 
fa^P) 



K -. / T/ \ r n i n 



(A8) 
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We next perform the integrals over the momenta. Using the integral representation of the S 
function, we have 
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(A9) 



where m = X^i=i m r We rewrite the full integral as 



n 



i i 

iVz! 2vri 



/ ^ n v A 



3/2 N 



/ K 



ra \ 
2tt/i 2 /3 J 



3/2 N 



(A10) 



x exp 



(3 ( 6 - B - E c0 - (/; (0) - Bi " ^ci) 
i=i 



Combining the factors resulting from the momentum integrals with those due to the fragment 
multiplicities, we define the translational Helmholtz free energies as 

3/2 N 



fi W) 
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( m N Ai \ ' 
\27rh 2 /3 J 



In (W) 



(All) 



where we recognize that mi = m^Ai and m = m^A , with the nucleon mass. We can 
then define the total Helmholtz free energy as 

k 

F n 09) = Ni if i 09) + /f a " S G») " A - S d ) - (/ trans 09) - Eco) , (A12) 

2=1 

where we have replaced the sum over the fragments with a sum over distinct fragments times 
their multiplicities, and write the density of final states as 



0O r . 



1 

2vri 



c+ioo 



dp exp [-0 {F n (/3) - E )\ 



(A13) 



where we have defined 

E = e -B . (A14) 

Up to this point, we have made several definitions but no approximations. This expression 
is exactly equivalent to that of Eq. ([9]). 
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